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15 Abstract. 

16 Results are presented on responses in 14-yr time series of stratospheric ozone and 

17 temperature from the Halogen Occultation Experiment (HALOE) of the Upper 

1 8 Atmosphere Research Satellite (UARS) to a solar cycle (SC-like) variation. The ozone 

19 time series are for ten, 20-degree wide, latitude bins from 45S to 45N and for thirteen 

20 “half-Umkehr” layers of about 2.5 km thickness and extending from 63 hPa to 0.7 hPa. 

21 The temperature time series analyses were restricted to pressure levels in the range of 2 

22 hPa to 0.7 hPa. Multiple linear regression (MLR) techniques were applied to each of the 

23 130 time series of zonally-averaged, sunrise plus sunset ozone points over that 

24 latitude/pressure domain. A simple, 1 1-yr periodic term and a linear trend term were 

25 added to the final MLR models after their seasonal and interannual terms had been 

26 determined. Where the amplitudes of the 1 1-yr terms were significant, they were in- 

27 phase with those of the more standard proxies for the solar uv-flux. The max minus min 

28 response for ozone is of order 2 to 3% from about 2 to 5 hPa and for the latitudes of 45S 

29 to 45N. There is also a significant max minus min response of order 1 K for temperature 

30 between 15S and 15N and from 2 to 0.7 hPa. The associated linear trends for ozone are 

3 1 near zero in the upper stratosphere. Negative ozone trends of 4 to 6%/decade were found 

32 at 10 to 20 hPa across the low to middle latitudes of both hemispheres. It is concluded 

33 that the analyzed responses from the HALOE data are of good quality and can be used to 

34 evaluate the responses of climate/chemistry models to a solar cycle forcing. 

35 
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35 1 . Introduction 

36 There have been a number of model studies focused on the response of middle 

37 atmosphere ozone to the 1 1-yr cycle of the solar ultraviolet (uv)-flux forcing (e.g., Garcia 

38 et al., 1984; Brasseur, 1993; Huang and Brasseur, 1993; Fleming et al., 1995; Harris el 

39 al., 1998; Shindell et al., 1999; Callis et al., 2000; Lee and Smith, 2003; Marsh et al., 

40 2003; Rozanov et al., 2004; Austin et al., 2007; Marsh et al., 2007; McCormack et al., 

41 2007; Nissen et al., 2007; Smith and Matthes, 2008; Austin et al., 2008). It is important 

42 to verily those modeled responses as a function of altitude, latitude, and season using 

43 high quality, global-scale observations of ozone. Analyses of the long-term satellite 

44 ozone data sets from the Stratospheric Aerosol and Gas Experiments (SAGE I/II) and 

45 from the series of Solar Backscatter Ultra-Violet (SBUV) instruments have yielded 

46 zonal-mean solar-cycle (SC) response patterns having maxima at middle latitudes from 

47 about 3 to 1.5 hPa (or ~40 to 45 km) [e.g., Wang et al., 1996; Chandra and McPeters, 

48 1994; McCormack and Hood, 1996; Lee and Smith, 2003; Soukharev and Hood, 2006; 

49 Randel and Wu, 2007]. However, the patterns for those responses do not agree so well 

50 with many of the simulations, which tend to peak at somewhat lower altitudes and lower 

51 latitudes. In addition, the SC amplitudes derived from the observations are about twice 

52 those from many of the models. On the other hand, Chandra and McPeters [1994] 

53 reported very good agreement between modeled and observed responses to the uv-flux 

54 variations associated with the 27-day solar rotation cycle. 

55 

56 Wang et al. [1996] found that the apparent SC response in the SAGE dataset was larger 

57 and occurred at higher latitudes for the northern hemisphere (NH) than for the southern 
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58 hemisphere (SH). Conversely, Lee and Smith [2003] found a SC response in the SBUV 

59 data that was larger in the SH, but they also cautioned that the observed responses may 

60 depend on the specific decade(s) for the given datasets. Those apparent differences have 

61 led researchers to consider the role of decadal-scale dynamical interactions for ozone, 

62 particularly for the middle and lower stratosphere [Nedoluha et al., 1998; Randel el al., 

63 2000; Tian et al., 2006; McCormack et al., 2007; Smith and Matthes, 2008]. Such 

64 considerations have prompted further model studies that include the effects of 

65 atmospheric interactions with the semi-annual oscillation (SAO) and quasi-biennial 

66 oscillation (QBO) cycles, as well as influences from volcanic events and from the El 

67 Nino/Southem Oscillation (ENSO). All of these forcings could alter a decadal-scale 

68 ozone response, perhaps even having a maximum that is in-phase with the SC for the 

69 stratosphere but located at middle latitudes rather than at tropical latitudes — or more like 

70 those obtained from the preceding data analysis studies [e.g., Matthes et al., 2004; Lee 

71 and Smith, 2003; Kodera and Kuroda, 2002], Hood [2004] summarized the nature of the 

72 model/data discrepancies for the 1 1-yr SC response. 

73 

74 A proper characterization of the SC response in ozone is also critical to a search for a 

75 slowdown for the decline of ozone during the 1 990s and/or for its recovery since then, in 

76 association with the effects of the mandated reductions for the release of 

77 chlorofluorocarbon compounds into the atmosphere [Reinsel, 2002; Reinsel et al., 2002; 

78 Newchurch et al., 2003; Steinbrecht et al., 2004a; Cunnold et al., 2004; Steinbrecht et al., 

79 2004b; Weatherhead et al, 1998; Randel and Wu, 2007], For example, Anderson and 

80 Russell [2004] and WMO [2007, Figure 1-12] report that the total chlorine derived from 
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81 measurements of HC1 obtained with the Halogen Occultation Experiment (HALOE) 

82 leveled-off after about 1997 and began to decline after 2001. Thus, it is reasonable that 

83 the trends for upper stratospheric ozone have been changing due to its chemical 

84 interactions, too [Siskind et al., 1998]. Rosenfield et al. [2005] reported that the chemical 

85 response of ozone ought to be asymmetric between the northern and southern 

86 hemispheres because of their differences for the distributions of the species responsible 

87 for ozone loss. 

88 

89 A proper characterization of the solar cycle or any other decadal-scale variation is also 

90 important for the determination of the ozone trend. For example, Soukharev and Hood 

91 [2006] analyzed for interannual, SC and trend effects, as part of their regression modeling 

92 of the SBUV, SAGE II, and HALOE ozone datasets. They focused primarily on the SC 

93 terms from their regression modeling, although they did show that the SC terms that they 

94 obtained from the SBUV dataset had response profiles at low latitudes that were similar 

95 for 1979-1991 and from 1992-2003 — decadal periods when the trends in upper 

96 stratosphere ozone were quite different. 

97 

98 Another important aspect of the fidelity of the observed SC response in ozone is how it 

99 relates to the SC response in temperature versus pressure, or T(p). Up until recently, that 

100 response in temperature had been obtained only with data from the NOAA operational 

101 sounding instruments, which provide relatively low vertical resolution radiance profiles 

102 for its retrieval in the stratosphere. Three distinct versions of the T(p) derived from the 

103 radiance data lead to SC responses that have a maximum at low latitudes, but which have 
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104 different magnitudes and are occurring over a different range of pressure-altitudes 

105 [McCormack and Hood, 1996; Scaife et al., 2000; Crooks and Gray, 2005]. There is also 

106 a problem with the drifting local time of the NOAA instruments, which could introduce 

107 long-term temperature variations because different times of the daily temperature cycle 

108 are sampled [Shine et al., 2008]. As a result, there is some uncertainty about which T(p) 

109 dataset is more appropriate for comparison with the observed SC ozone responses. 

110 

111 This paper contains results of analyses for the decadal-scale responses and trends in 

1 12 ozone and T(p) based on data for 1991-2005 from the single, well-calibrated HALOE 

113 instrument. In particular, it is focused on the characterization of the response of its upper 

114 stratospheric ozone and temperature to the 1 1-yr SC, keeping in mind that the HALOE 

115 measurements span only one complete SC. For this reason the results herein may be 

1 1 6 considered somewhat exploratory, rather than definitive. Because the trends in upper 

117 stratospheric ozone due to its interaction with chlorine were small during the time span of 

118 the HALOE measurements [Rosenfield et al, 2005; Yang et al., 2006], it was felt that one 

119 could have better success with isolating a SC-like response for the time period of 

120 HALOE rather than for the datasets that included the 1 980s, too. Also different from the 

121 two previous solar maxima, there was no interference from a volcanic eruption near the 

122 solar maximum of 2001 [Solomon et al, 1996]. Because the calculated percentage SC- 

123 like response in ozone is much larger than its percentage response in T(p) it ought to be 

124 easier to isolate an SC-like signal in ozone. Nevertheless, one must still be aware that it 

125 is easy to confound the effects of a SC term and a trend term, especially when adjacent 
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126 points in the time series of the zonal means are correlated, and/or when the end portions 

127 of the time series are anomalous (e.g., Tiao et al. [1990]). 

128 

129 Section 2 reviews the attributes of the HALOE ozone profiles, along with the time series 

130 approach to their analysis. Section 3 includes the findings from those analyses and 

131 describes the important similarities and differences with the SC responses from the 

132 several other published, observed and modeled studies. Section 3 also relates the SC-like 

133 responses from analyses of the HALOE temperature time series [Remsberg and Deaver, 

134 2005; Remsberg, 2007; Remsberg, 2008] and reviews how it compares with results from 

135 zonal-mean models and observations. In general, there is good agreement with the model 

136 predictions of the SC responses in ozone and T(p) of Brasseur [1993], Marsh et al. 

137 [2007], and Austin et al. [2008] versus those found in the HALOE data, except in the 

138 tropical upper stratosphere. Linear trends from the ozone analyses are presented in 

139 Section 4, along with similar findings for the temperature from the more restricted 

140 pressure-altitude range of 2 to 0.7 hPa. Section 4 also comments on a tropical to 

141 subtropical anomaly in the SC-like responses of the HALOE ozone in the lower 

142 stratosphere, at least as compared to the modeled responses of ozone to just the 1 1 -yr uv- 

143 flux forcing. Summary findings are in Section 5. 

144 

145 

146 
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146 2. Approach 

147 a. Data characteristics 

148 HALOE is a solar occultation instrument that operated successfully on the Upper 

149 Atmosphere Research Satellite (UARS) from October 1991 through late November 2005. 

150 A description of the HALOE experiment is given in Russell et al. [1993], and a first 

151 characterization of its ozone profiles is provided in Briihl et al. [1996]. An update of the 

152 quality of the current operational ozone dataset is given in Randall et al. [2003] for the 

153 so-called Version 19 (or VI 9) Level 2 profiles that are used for the present study. The 

154 HALOE Project Team monitored the performance of the HALOE instrument over its 

155 mission lifetime, and they found no long-term degradation of its measurement 

156 characteristics that would impact the determination of a SC or long-term trend in ozone 

157 [Gordley et al., 2006; Hervig et al., 2007]. 

158 

159 The HALOE ozone profiles were obtained with excellent signal-to-noise (S/N) and with 

160 no a priori constraints for their retrievals. They have a vertical resolution of about 2.3 

161 km and are registered in pressure-altitude, in particular as derived from a retrieval of its 

1 62 own temperature profiles from about 4 hPa (near 38 km) to 0.004 hPa (near 85 km). The 

163 vertical resolution and the measurement and retrieval sensitivity of the HALOE ozone 

164 ought to be providing an accurate SC response profile, as compared with that of the more 

1 65 standard SBUV sounding technique, for example. 

166 

167 The T(p) profile information in the HALOE dataset below about the 4-hPa level was 

1 68 obtained from the daily 12Z operational analyses provided to the UARS Project by the 
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169 NOAA/Climate Prediction Center (CPC) and has relatively low vertical resolution. The 

170 HALOE T(p) information above that level has a vertical resolution of 3-4 km and is 

171 representative of the local times of the occultation measurements, both of which also 

172 provide advantages when compared with the CPC T(p) data [Remsberg, 2007; Remsberg, 

173 2008]. Thus, the HALOE ozone and T(p) datasets are compatible above the 4-hPa level 

174 for obtaining SC-like responses from the upper stratosphere to the lower mesosphere and 

175 for relating them to each other. 

176 

177 The HALOE measurements are for sunrise (SR) and sunset (SS), and there are significant 

178 variations in ozone at those two times in the low to middle mesosphere due to rapidly- 

179 changing, photochemical effects. First-order, diurnal corrections are included in the 

1 80 forward model for HALOE VI 9 ozone above about 37 km to account for those changes 

181 along the limb tangent paths for the retrieved profiles. Natarajan el al. [2005] re- 

1 82 evaluated the adequacy of those corrections and developed further modifications that 

1 83 ought to be applied to updated HALOE algorithms for the retrieval of its mesospheric 

184 ozone, especially for SR conditions and above about 61 km (or about 0.2 hPa). In 

1 85 addition, Marsh et al. [2003] showed that seasonal cycles for upper mesospheric ozone 

1 86 are small at SR, but not at SS, and that it is not a good idea to combine their SR and SS 

1 87 data for the mid mesosphere when searching for its SC response. For these reasons the 

188 analyses for a SC response in this study are restricted to HALOE VI 9 profile segments 

1 89 below the 0.7-hPa level (or about 53 km). 

190 
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191 Figure 1 shows the measurement pattern of the tangent point locations from HALOE for 

192 each day of 2001. One can see that the measurements were less frequent poleward of 

193 about 50 degrees latitude and are missing at 40 and 50 degrees during summer. Prior to 

194 1996, the summertime and mid latitude gaps are not as large because of the greater power 

195 availability for instrument operations during the early years of the UARS mission. The 

196 intervals for adjacent samples within a latitude zone from solar occultation measurements 

197 vary from a few days to more than a month. Although rather infrequent and variable, 

198 such sampling intervals are still adequate for resolving the seasonal and longer-term 

199 variations from the total data record and with sufficient accuracy. 

200 

201 b. Analysis procedure 

202 Some initial studies of the seasonal and longer-period variations in HALOE ozone were 

203 conducted for the lower stratosphere by Remsberg et al. [2001] using multiple linear 

204 regression (MLR) techniques, and they found that annual and quasi-biennial (QBO) terms 

205 accounted for almost all of its variations. In the current study MLR was applied to 

206 HALOE time series of averages of the ozone points for a given pressure layer and for 20- 

207 degree wide latitude bins centered every 10 degrees from 45N to 45S. Latitude bins of 

208 this width provide a high probability of having enough samples to adequately represent 

209 the zonal mean. The MLR technique also accounts for the fact that the data points from 

210 occultation measurements are inherently non-orthogonal. The HALOE Level 2 profiles 

211 are tabulated with a spacing of 300 m. However, the precision of the data points for an 

212 ozone time series was improved by linearly integrating the individual ozone profiles in 

213 log pressure over “half-Umkehr layers” (about 2.5 km thick or very near to the vertical 


10 



214 resolution of the retrieved profiles). The resulting ozone amounts for those layers are 

215 given in terms of Dobson Units (DU). Table 1 defines the half-Umkehr layers for this 

216 study. The total pressure-altitude range is 63.3 hPa to 0.70 hPa or a range of 13 half- 

217 Umkehr layers labeled from 4L to 10L, respectively. 

218 

219 Other researchers have reported on how the SC ozone response in the lower stratosphere 

220 is confounded with and possibly dominated by decadal-scale dynamical effects [e.g., 

221 Dunkerton and Baldwin, 1992; Salby et al., 1997]. Thus, the current study is focused on 

222 the middle and upper stratosphere, where the dynamical effects are weaker and hopefully 

223 more easily resolved. In most other respects the analysis approaches in Remsberg et al. 

224 [2001] and especially in Remsberg [2007; 2008] were followed. Annual (AO) and semi- 

225 annual (SAO) cycles were fit to the time series, and the time series residuals were 

226 analyzed for their remaining interannual structure. Weak quasi-biennial or QBO-like (28 

227 mo.) and sub-biennial or IA (21 mo.) terms were also found consistently in the upper 

228 stratosphere. The QBO periods became more variable in the middle and lower 

229 stratosphere (25 to 3 1 months), in accord with the power spectral analysis results of 

230 HALOE ozone by Witte et al. [2008]. The corresponding subbiennial periods changed, 

23 1 as well. An overall objective was to account for any significant, periodic structure in the 

232 time series before analyzing for the underlying SC-like and trend terms. 

233 

234 The average separation in time for the points in the HALOE time series is about 25 days. 

235 Even so, because zonal-average ozone values at time n have considerable memory of the 

236 atmospheric state at time n-1, the time series have autoregressive characteristics. The 
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237 appropriate MLR result for that situation is obtained by a two-step process (see Appendix 

238 A of Tiao et al. [1990]). Initially, a model of the form, 

239 

240 03 (n) = a + bX + e sin (Z) + f cos (Z) + ... + N(n), (1) 

241 

242 was used, where 03 (n) is the sum of the components of ozone at the n th point in the time 

243 series. The model has a constant term ‘a’ and a linear term in X = (t n - ti) / T, where t n is 

244 the time of the n th observation point, ti is the first point in the time series, T is the total 

245 length of the time series, and b is the coefficient of the linear term. The model has 

246 seasonal and longer-period terms that are represented by Fourier components in Z, where 

247 Z = 27rt n / P, and P is the period of a given cycle. Periodic terms that were considered are 

248 AO, SAO, QBO-like, sub-biennial (or IA), and an 1 1-yr (4017-dy) term. Note that the 

249 seasonal terms were almost always found to have a probability of at least 99% of being 

250 present. The QBO, IA, 1 1-yr, and linear terms were allowed to have lower probabilities, 

25 1 in order to accommodate their existence in the time series across all latitudes and 

252 pressure-altitudes. Finally, the term N(n) is the autoregressive (AR) noise residual that is 

253 given by 

254 

255 N(n) = cp N(n-l) + E(n), (2) 

256 

257 for a first order (or AR1) process. The factor cp is the autocorrelation of the noise residual 

258 and E(n) is the white noise component. Thus, the MLR method consists of a fit for all the 

259 terms of Eq. (1) and the generation of its model residuals N(n). The autocorrelation 
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260 coefficient cp for this model is found from the residuals. The model terms of Eq. (1 ) are 

261 then transformed to those of Eq. (3), 

262 

263 03 (n) - tp 03(n-l) = a [1 - <p] + b [X(n) - cp X(n-l)] + (3) 

264 e [sin(Z(n)) - cp sin(Z(n- 1 )] + f [cos(Z(n)) - (p cos(Z(n- 1 ))]+.. .+ E(n), 

265 

266 and the data time series is refit to get new coefficients. Initially, (p was set to zero for the 

267 MLR analyses using the relevant terms of Equation (1), cp was calculated from the model 

268 residuals, and then final coefficients for the terms were obtained with Equation (3). 

269 Because MLR provides for a combined fitting of all the model terms, the absolute ozone 

270 uncertainties are nearly the same for each term of the model. Percentage deviations and 

27 1 the probabilities of their respective terms depend on their amplitude. 

272 

273 Figure 2 is an example of a HALOE data time series (Layer 8U centered at about 42 km 

274 and at the latitude of 25N), where the SR and SS data have been fit with a regression 

275 model (the solid oscillatory curve) that contains the seasonal, AO and SAO terms, a 

276 QBO-like (853-dy) term, a subbiennial (640-dy) term, and an 1 1-yr term. The horizontal 

277 straight line is the constant term from the model, and all terms are highly significant. An 

278 1 1 -yr term was fit to the time series rather than regressing against a proxy term for the 

279 uv-flux, to determine whether there might be other decadal-scale terms that were 

280 affecting the ozone. The 1 1-yr response in Figure 2 is essentially in-phase with the time 

281 series of the solar flux for SC 22 and 23 (see following subsection, too). In other words, 

282 there is a maximum in the ozone in this case that is very near to the time of a maximum 
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283 in the solar flux (i.e, 0.9 yr from January 1991 and 2002), as expected from increases for 

284 the uv-photodissociation of molecular oxygen [Brasseur and Solomon, 1984], One can 

285 also see that the maxima for the annual cycles do not exceed the horizontal line in 1 996 

286 and 1998, or near solar minimum; its 1 1-yr, maximum minus minimum response is 2.7%. 

287 For this layer the amplitudes of the seasonal and interannual terms are small and of the 

288 order of the amplitude of the estimated 1 1 -yr or SC-like term. At this point it is noted 

289 that the zonal mean distribution of ozone can have significant meridional variations at 

290 mid and higher latitudes of the middle and upper stratosphere, yet Eq. (3) does not 

291 contain a term to account for variance with latitude. Such variations are embedded 

292 mainly in the coefficients of the seasonal terms. Since the latitude bins are relatively 

293 narrow (20 degrees wide) compared with the meridional gradients in the zonal-average 

294 ozone, the data time series provide good estimates of the true seasonal coefficients for 

295 zonal-mean ozone. The aperiodic ENSO or volcanic forcings were not modeled; their 

296 effects may be important for ozone in the lower stratosphere (see later). 

297 

298 The model fit in Figure 2 does not include a linear term, although that term is included 

299 for the MLR results of Sections 3 and 4. Figure 3 is a plot of the residuals for the ozone 

300 model of Figure 2. It shows that there is no remaining apparent periodic structure, which 

301 is an important acceptance test for the set of terms of the final MLR model. There are 

302 also no significant, non-random autocorrelations in the residuals for this final model (not 

303 shown, but see Remsberg et al. [2001] for a diagnostic plot of this statistical test). The 

304 solid horizontal line is the linear fit to the residuals; and in this case it is not significantly 

305 different from zero. A polynomial trend term may account for a chemical response of 
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306 ozone to the long term changes in total reactive chlorine [e.g., Randel and Wu, 2007; 

307 Yang et al., 2006], but it was not included in the MLR models because reactive chlorine 

308 is not the dominant loss mechanism for ozone throughout the stratosphere. Furthermore, 

309 the time series of the residuals in Figure 3 does not indicate the presence of such an 

310 underlying, polynomial trend term due to that chemical loss. However, since total 

311 reactive chlorine reached its peak value in the upper stratosphere in the late 1 990s [ WMO ; 

312 2007], any enhanced chemical loss at that time would tend to reinforce the minimum 

313 response of ozone to the SC forcing. Thus, the 2.7% max minus min, SC-like response of 

314 the MLR model fit in Figure 2 may overestimate the amplitude of the true SC effect on 

315 ozone. 

316 

317 Ozone is in photochemical equilibrium in the upper stratosphere, and its seasonal 

3 1 8 variations are anti-correlated with the corresponding seasonal variations for HALOE 

319 T(p), as shown in Figure 4 for 25N and 2 hPa. A linear trend term is included in its 

320 model, and its diagnosed value is -1.1 K/decade. The 1 1-yr variation in T(p) is much less 

321 apparent (max minus min is 0.6 K). That percentage response (0.23%) is smaller than for 

322 the ozone but in accord with at least some model predictions [Austin et al., 2008], Even 

323 so, the diagnosed phase for the maximum of the 1 1-yr term in T(p) occurs 1.7 years prior 

324 to solar max and leads the 1 1-yr response of ozone by 2.6 years based on its MLR model 

325 in Figure 2. The adjusted, in-phase, max minus min SC response is more like 0.3 K for 

326 T(p). Therefore, it may be that there are other effects contributing to the analyzed, 1 1-yr 

327 response in T(p) for this layer and latitude (see later in Section 4). 

328 
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329 Finally, one should note that the SR and SS samples in the HALOE time series of Figures 

330 2 and 4 tend to alternate, due to the sequential crossings of the orbital sampling patterns 

33 1 for the SR and SS measurements of Figure 1 . This character for the sampling leads to 

332 largely real, diurnal variations in both the ozone and temperature data in the mesosphere. 

333 Its occurrence enhances the ‘noisy’ appearance of the HALOE time series residuals after 

334 removal of the seasonal and interannual terms. In addition, those alternating SR/SS 

335 residuals are negatively correlated, obscuring the otherwise positive autocorrelation of 

336 the residuals at lag-1 and higher lags. A first-order, average adjustment for those small, 

337 SR/SS differences was obtained for these higher altitude data; mean values of the 

338 separate time series of the SR and of the SS data were determined and half the difference 

339 of those two means was then applied to the respective SR and SS points. This step 

340 reduces the short-term, point-to-point, fluctuations in the time series and increases the 

341 significance of the AR1 term and, thus, the other small amplitude terms for the regression 

342 models of Figures 2 and 4, but it has essentially no effect on the phases of the periodic 

343 terms. Therefore, for each time series the adjusted SR and SS points were combined, 

344 giving about 200 points in all. A clearer example of this effect for the temperature tides 

345 of the mesosphere can be found in Remsberg [2008, their Figures 1 and 2]). On the other 

346 hand, almost no systematic, SR/SS differences were found for the time series of the 

347 HALOE ozone in the mid and lower stratosphere, as expected based on model 

348 calculations of its weak diurnal variations. 

349 

350 
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350 c. Signatures of interannual terms from the MLR models 

351 Traditional proxies for the QBO term were tried during the early development of the 

352 MLR models for this study. However, QBO cycles based on the tropical lower 

353 stratospheric winds did not provide for a highly significant fit of the interannual 

354 variations in upper stratospheric ozone, particularly for the subtropics of the NH. Similar 

355 deficiencies for those QBO-proxies were reported for time series analyses of SC effects 

356 in temperature by Crooks and Gray [2005]. In the present study the QBO-like terms that 

357 were fitted to the ozone residuals for the low latitudes had significant changes of phase 

358 from 50 to 10 hPa, not unlike that reported by Pascoe et al. [2005] for the zonal winds. 

359 Thus, a QBO wind proxy that is averaged over that pressure range may not be so 

360 representative of the true forcings of that term. Furthermore, the effect of a QBO forcing 

361 in ozone is a bit more complicated because it also depends on the zonal -mean ozone 

362 distribution and its gradients, which change in magnitude and sign for the vertical range 

363 of the subtropical stratosphere. Instead, the periods and phases of the dominant 

364 interannual terms were obtained for the present study from a Fourier analysis of the 

365 points of each of the time series. 

366 

367 Figure 5 shows the time series of HALOE ozone observations for 25N and Layer 6L 

368 (1 1.2 - 15.8 hPa). The combination of the AO and SAO terms of the MLR model 

369 provides a very good fit to the seasonal variations for this layer. Amplitudes for the 

370 interannual and SC-like terms are weaker, although it is obvious that there are significant 

371 interannual variations in the data. Those variations are fit well by a combination of a 26- 

372 month QBO-like term and a 22-month subbiennial term. 
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373 


374 The ozone of Figure 5 for 1991-92 is clearly elevated compared with that of the other 

375 years, most likely a result of physical processes associated with the major volcanic 

376 aerosol layer from the Mt. Pinatubo eruption of June 1991. In fact, the ozone is 

377 somewhat anomalous during this period for other low to middle stratospheric layers in 

378 the tropics and subtropics of both hemispheres, and Section 4 contains some discussion 

379 about possible causes of the excess ozone of 1 991-92. A model fit of the entire time 

380 series of Figure 5 was attempted, including both a linear trend term and a periodic term of 

381 1 1-yr period. Ftowever, the excess ozone of 1992 represents an “endpoint anomaly” for 

382 the analysis. It leads to a rather large, decreasing trend and affects the nature of the 

383 associated 1 1-yr term. Instead, for those layers that showed anomalies in this early 

384 period, it was found that much more reasonable results were obtained for those two terms 

385 by starting the MLR model in early September 1 992 or in January 1 993. That 

386 modification was applied to the analyses of the layers 6U through 4L. The solid 

387 oscillating curve in Figure 5 is the MLR model fit to the data using all terms, but 

388 beginning in September 1992. Note also that the solid curve actually begins several 

389 points later because of the need to allow for the autoregressive nature of the time series. 

390 

391 The 1 1 -yr term of the model fit in Figure 5 has a “max minus min” value of 2.7%, but the 

392 phase of its maximum occurs 1.0 year prior to January 1991 (or 2002). A decreasing 

393 trend of 5. 1 %/decade was also realized from the analysis. The results for the 1 1 -yr terms 

394 from each of the 130 ozone time series are presented in Section 3. Results from similar 

395 analyses of the HALOE T(p) are given there, as well, but only for the pressure layers of 


18 



396 0.7 to 2.0 hPa. Findings for the linear trend terms are presented and discussed in Section 

397 4. The rest of this subsection considers the QBO-like terms and their interactions with 

398 the seasonal terms. 

399 

400 A typical example of the zonal mean ozone distribution is shown in Figure 6 from a 

401 sequence of HALOE SS observations that eventually covered most latitudes during the 

402 period of February 26 to April 5, 1995. The QBO cycle has its largest effect on the zonal 

403 mean ozone at the tropical latitudes and at those pressure altitudes where the vertical 

404 ozone gradient is large. In fact, one can check on the realism of the current QBO signal 

405 in ozone by examining its amplitudes from the analyses of the 130 separate time series, as 

406 shown in Figure 7. QBO amplitudes are nearly hemispherically-symmetric throughout 

407 the stratosphere. The largest amplitudes (in %) occur in the tropical stratosphere near 7 

408 hPa and 30 to 40 hPa. Photochemistry controls the distribution of ozone above about the 

409 5-hPa level and thus dampens the effect of the QBO on the vertical ozone gradient of the 

410 upper stratosphere. Its amplitudes are also relatively weak at about 20 hPa for the 

411 latitudes of 15 to 45 degrees of both hemispheres, which is where the meridional 

412 gradients of ozone are weak, as well (c.f., Figures 6 and 7). 

413 

414 Dunkerlon [2001] explained that a subbiennial term occurs in ozone due to the difference 

415 interaction of the annual and QBO cycles. Figure 8 is a plot of the amplitudes of the 

416 subbiennial terms. Again there is good continuity for the features of their distribution, 

417 and their amplitudes decrease with altitude in the region of transition from dynamical to 

418 photochemical control. Maximum amplitudes occur in the subtropics of the middle 
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419 stratosphere, where the combined effects due to the annual cycle and the QBO are largest 

420 for ozone. In the tropics near 30 hPa and 7 hPa the subbiennial amplitudes are weaker 

421 than that for the QBO forcing alone because the annual cycle is weak at low latitudes. 

422 The subbiennial amplitudes are somewhat larger in the northern than the southern 

423 hemisphere subtropics at about 10 hPa. The annual cycle amplitude is also larger for the 

424 northern hemisphere in this region, presumably an indication of the zonal mean effects of 

425 the more vigorous and sustained net transport during winter and springtime for the 

426 northern hemisphere. The region of maximum values near 50 hPa and 30 degrees of 

427 latitude is also where the annual cycle in ozone is relatively large. 

428 


20 



428 3. Solar Cycle Responses in HALOE Data 

429 a. Ozone 

430 Initially, the time series of F10.7 cm solar radio flux data was employed as a proxy for 

43 1 the SC term, but its fit to the HALOE ozone was not significant at the 90% confidence 

432 level. That proxy term exhibited shorter period structure that was not present in the 

433 ozone time series of HALOE even after applying a smoothing to the daily values of the 

434 proxy with an 81-dy running mean to avoid the effects of the 27-dy solar rotation cycle in 

435 the daily fluxes. As mentioned earlier, there may also have been a confounding of the 

436 effects of the solar uv-flux with those from reactive chlorine during 1 991-2005, at least in 

437 the upper stratosphere. Therefore, it was decided that it would to be instructive to fit an 

438 1 1-yr periodic term to the time series but to allow its phase to be determined simply from 

439 its best fit to the data. 

440 

441 Figure 9 shows the max minus min ozone values (in percent) for its 1 1 -yr terms over the 

442 domain of latitude and pressure-altitude. Maximum responses of order 2 to 3% were 

443 found in the upper stratosphere at subtropical to middle latitudes, and they are somewhat 

444 asymmetric between the two hemispheres. The dark shading shows where the 1 1-yr 

445 terms have a confidence interval (Cl) of greater than 90% for their partial rank order 

446 correlations [Remsberg et al., 2001], and the lighter shading indicates where they have a 

447 Cl value between 70 and 90%. Max minus min values of up to 2% extend through the 

448 tropics near 3 hPa, and they remain significant at the 90% level. This finding indicates 

449 that the sampling from HALOE was adequate for representing the zonal average ozone in 

450 the tropics. In the lower stratosphere the largest percentage responses occur in the 
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45 1 subtropics of the southern hemisphere and in the tropical latitudes (although larger than 

452 expected at 50 hPa). Smaller and much less significant responses occur in the tropics 

453 near 1 hPa (0.5%) and from 7 to 30 hPa (1%), at low to middle latitudes of the southern 

454 hemisphere near 7 hPa (1 to 1.5%), and at middle latitudes of the northern hemisphere at 

455 40 hPa (1%). However, it is stressed that the estimates of significance assume that all the 

456 relevant structure in the time series has been accounted for. 

457 

458 To see how well the diagnosed, 1 1-yr terms are related to the solar cycle flux, their 

459 phases are given in Figure 10 in units of years from January, 1991. Due to the fact that 

460 the maximum for SC 22 was broad and nearly flat over 2 years, it was assumed that the 

461 1 1-yr signal in ozone was in-phase if it was within ±1 yr of solar uv-flux maximum (the 

462 regions delineated by shading in Figure 10). SC 22 had its final peak in early 1 992 

463 followed by a sharp decline, and it has been argued by some that phase maximum for the 

464 uv-flux occurred at that time rather than a year earlier in January 1991. But in order to be 

465 consistent with the estimates of the phase of the SC maximum for T(p) of Rems berg 

466 [2007; 2008], that January 1991 time of reference was kept for the ozone results, too. 

467 Generally, the 1 1-yr terms for ozone are in-phase with that of the uv-flux, with the 

468 exception of the upper stratosphere at tropical latitudes and near 30 hPa from the tropics 

469 to northern hemisphere middle latitudes. It is presumed that the lower stratosphere 

470 anomaly is indicative of an interannual, dynamical forcing mechanism that is dominating 

471 the effects of the uv-forcing. In fact, Marsh and Garcia [2007] reported ozone at 52 hPa 

472 that was anti-correlated with a lagged ENSO index; they found that this dynamical 

473 forcing tended to reduce the response of ozone to a solar forcing in the lower stratosphere 
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474 during 1991 to 2005. Still, the results from Figure 10 indicate that the presumed ENSO 

475 forcing does not seem to be affecting southern hemisphere ozone in the same way. 

476 

477 Figure 10 is showing that the 1 1-yr term has a maximum that lags that of the uv-flux by 

478 several years in the tropical stratosphere from 4 to 1 hPa. This anomaly may be an 

479 indication of an enhanced chemical loss of ozone due to reactive chlorine, sometime after 

480 solar minimum. It is noted that such losses for annual-averaged ozone ought to be most 

481 clearly evident in the upper tropical stratosphere because that is where the effects of the 

482 chlorine would began to compete with those of the uv-flux [Brasseur and Solomon, 

483 1984], In support of this idea, it is noted that the largest values of total chlorine, as 

484 inferred from the HALOE measurements, occurred from 1997 to 2002 (see Figure 1-12 

485 of WMO [2007]). In this instance the inclusion of a simple linear, rather than a 

486 polynomial, trend term in the MLR models may not adequately represent the changing 

487 effects of the chlorine. Indeed, if this process is contributing to the out-of-phase, 1 1 -yr 

488 ozone response, it would be unique to the time span of the HALOE observations; that is, 

489 total chlorine reached its maximum values in the upper stratosphere in the middle of its 

490 14-yr period. 

491 

492 b. Temperature 

493 A test of the fidelity of the ozone results is an assessment of the analogous 1 1 -yr response 

494 in temperature from HALOE, even though model simulations indicate that the percentage 

495 change for T(p) ought to be much smaller than that for ozone (e.g., Smith and Matthes 

496 [2008]). Remsberg and Deaver [2005] analyzed for that response using 10-degree 
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497 latitude bins for the upper stratosphere, as well as the mesosphere. Their MLR analysis 

498 approach was essentially the same as that used here. However, they were not confident 

499 of their results at and below the 3-hPa level because of concerns about the long-term 

500 accuracy of the temperature time series from the NOAA/CPC to which the HALOE T(p) 

501 retrievals were merged. As a result, the HALOE findings for the 1 1 -yr and trend terms in 

502 T(p) are only appropriate above the 3-hPa level for comparisons with those of ozone 

503 [Remsberg, 2007; 2008]. 

504 

505 The 1 1-yr temperature responses were obtained for 20-degree wide latitude bins in the 

506 present study, in order to be compatible with the ozone results. Those responses are 

507 shown in Figure 1 1 along with their phase deviations (in years) from January 1991 in 

508 Figure 12, but only for the pressure levels of 0.7 to 2.0 hPa. The 1 1 -yr temperature 

509 maximum is in phase (within ±1 yr) with the uv-flux forcing from about 25S through 

510 15N. There are transitions to somewhat out-of-phase values at 35S and at 25N, but 

511 where the responses in Figure 1 1 are weaker and not very significant. 

512 

5 1 3 The HALOE 1 1 -yr response is smaller than reported from several analyses of the long- 

514 term temperature data sets from the operational satellites [Hood, 2004], On the other 

5 1 5 hand, the magnitude and variation with latitude and altitude for the FLALOE terms at 2 

516 and 1 hPa agree very well with the observed responses based on the recent findings of 

517 Crooks and Gray [2005, their Fig. 2] for the period 1979-2001, as based on the ERA-40 

5 1 8 reanalysis dataset. 
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519 The response in HALOE T(p) near the tropical stratopause is about two-thirds that found 

520 from the simulations of Brasseur [1993, their Figure 14b] and about half that from the 

521 fixed dynamical heating model of McCormack and Hood [1996]. However, the variation 

522 with latitude of the HALOE response compares favorably with the model results of 

523 Marsh et al. [2007] and of Smith and Matthes [2008], at least after a scaling of the latter 

524 results to the actual max minus min values of the solar flux. The solar cycle response 

525 from the AMTRAC model of Austin et al. [2008] most closely matches that reported 

526 here. 

527 

528 The temperature response in the HALOE data is in-phase with the solar uv-flux near the 

529 tropical stratopause. Thus, it is the HALOE ozone response that appears anomalous in 

530 that region. However, the amplitude of the 1 1-yr term in ozone is becoming small and is 

531 not significant at the Equator at 1 hPa. It is also noted that there was a peak in the uv- 

532 flux in about January 1 992, rather than January 1991. An allowance for that fact reduces 

533 the phase discrepancy for ozone by 1 year without affecting the phase agreement of the 

534 temperature response noticeably. 

535 

536 c. Solar cycle (or SC-like) responses of HALOE ozone and temperature 

537 In order to obtain estimates of the actual response of atmospheric ozone and temperature 

538 to a proxy for the direct forcing from the uv-flux, their sets of responses from Figures 9 

539 and 1 1 have been adjusted by the amount that they were not exactly in-phase from 

540 Figures 10 and 12. Those adjustments were made by multiplying the responses by the 

541 factor, cos [2n p / 11], where p is the deviation of phase of the maximum (in years) from 
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542 January 1991. Even after applying the adjustments, there is good continuity among the 

543 ozone responses from the 130 separate time series. 

544 

545 The separate, 1 1-yr responses for ozone of Figure 9 have been averaged for the latitude 

546 bins that extend from 25 S and 25N, and those profiles are shown in Figure 13. The 

547 corresponding, SC-like response profile is shown as well, but there is very little 

548 difference between the two. There is good similarity between the results in Figure 13 and 

549 the SC response profile obtained from the seasonally-averaged HALOE data points over 

550 the same latitude range in Soukharev and Hood [2006, their Figure 14], if one excludes 

551 their analyzed value at about 35 km. The magnitude and variation of the HALOE 

552 response profiles agree very well with the response at low latitudes from the 

553 representative, zonal-mean models of Brasseur [1993, his Figures 9a and Figure 14a] and 

554 of Huang and Brasseur [1993, their Figure 8a]. They also agree qualitatively with the 

555 results of Austin et al. [2007] and of several other models, as shown in Soukharev and 

556 Hood [2006, their Figure 14]. In addition, there is no discrepancy for this HALOE ozone 

557 response across most of the stratosphere versus the responses from the chemistry/climate 

558 models, AMTRAC and WACCM, in Austin et al. [2008]. 

559 

560 Lee and Smith [2003] analyzed both the SAGE II and SBUV data for a solar cycle 

561 response. Their results from SAGE II are similar to the zonal mean contour values of 

562 Figure 9, at least for the upper stratosphere at most latitudes; their results from SBUV do 

563 not agree as well though, especially in the tropics (see also the analyzed results from 

564 SAGE and SBUV in Soukharev and Hood [2006]). Lee and Smith [2003] conducted 
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565 model studies of the SC response, but modified to account for the relative phases of the 

566 QBO over specific decades. The interaction with the QBO in their model leads to larger 

567 positive responses in the subtropical upper stratosphere of both hemispheres, but their 

568 results are also sensitive to the QBO indexes that they employed. However, it is noted 

569 that for the present HALOE analyses the interactions with the QBO and its associated 

570 subbiennial term have been accounted for to first order based on the structure of those 

571 terms in the HALOE time series. Austin et al. [2007] were able to simulate the minimum 

572 ozone response at about 20 to 30 hPa, but not the analyzed, sharply increasing response at 

573 about 50 hPa in Figure 13. The model simulations reported in Marsh et al. [2007, their 

574 Figure 8] indicated a positive ozone response near 50 hPa at low latitudes for the period 

575 of 1 979 to 2003, but not for the longer period of 1 950 to 2003. It is presumed that the 

576 large response diagnosed from the HALOE data is not strictly related to the solar cycle. 

577 

578 The SC-like, max minus min responses for ozone (in percent) and for temperature (in K) 

579 are summarized in Table 2 for the region of the stratopause and for three distinct latitude 

580 zones, 45S to 35S, 15S to 5N, and 25N to 45N. For the tropical latitude zone the adjusted 

581 ozone response is small and not significant, while the associated temperature response is 

582 1 .0 K and highly significant. For the middle latitude zones the ozone response varies 

583 between 1 .9 and 2.8% and is highly significant, while the temperature response decreases 

584 to near zero. 

585 
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585 4. Trends in HALOE ozone 

586 Figure 14 shows the linear trends (in %/decade) obtained from the MLR analyses of the 

587 HALOE ozone. Darker shading is where those terms have a Cl of greater than 90%; 

588 lighter shading denotes a Cl between 70 and 90%. The trends are near zero in the upper 

589 stratosphere across most of the latitude domain and for this 14-yr time span, when the 

590 total chlorine levels had leveled off and started to decline slowly [ WMO , 2007]. The 

591 trends are decreasing by 4 to 6 %/decade across most latitudes of the middle stratosphere, 

592 and they are significant. There are net increases in the ozone of the lower stratosphere in 

593 the subtropics of both hemispheres, but slight decreases in the tropics. Those trends are 

594 not significant, however. Overall, there is very good continuity in the patterns of the 

595 trends for the latitude and pressure-altitude domain of Figure 14. 

596 

597 The trends shown in Figure 14 appear to be somewhat weaker than those reported from 

598 F1ALOE, SAGE II, and lidar in WMO [2007]. However, those analyses were for data 

599 time series on constant altitude surfaces, and it is noted that the effects of a long-term 

600 cooling of the atmosphere will lead to a lowering in altitude of a pressure surface and its 

601 ozone [Rosenfield el al., 2005]. For example, the temperature trends reported by 

602 Remsberg [2007; 2008] for 25S to 25N are of order -0.8 to -1.2 K/decade for the pressure 

603 surfaces of 0.7 to 2.0 hPa, respectively, or for those levels where it is possible to conduct 

604 reliable analyses for a temperature trend from retrievals of the HALOE measurements 

605 themselves. Thus, a negative ozone trend at an altitude level will appear as a weaker 

606 trend at a corresponding pressure level. 

607 


28 



608 There may also be slight biases in the analyzed trends from the HALOE zonal-averaged 

609 ozone versus pressure profiles at all latitudes for the layers 7L to 4U. This prospect is 

610 because the HALOE T(p) profiles for those layers were based on the data from the 

611 NOAA operational temperature sounders that were not corrected for the trends in the 

612 atmospheric CO 2 , as they affect the vertical weighting functions of their several 

613 measurement channels [Shine et al., 2008]. However, any such effects should not alter 

614 the trends with latitude in Figure 14 for a given pressure layer. 

615 

616 A number of investigators have wondered to what extent the eruption of Mt. Pinatubo 

617 may have affected the SC-like and trend terms for ozone [e.g., Lee and Smith, 2003; Al- 

618 Saadi et al., 2001]. Furthermore, Crooks and Gray [2005] point out that there were 

619 interactions between the QBO and ENSO that were altering the zonal winds of the 

620 subtropical middle stratosphere both prior to and just following the eruption of Mt. 

621 Pinatubo. Recall from Figure 5 that there was an excess of ozone in layer 6L at 25N 

622 during late 1991 to mid 1992. That excess may be related to atmospheric effects due to 

623 the presence of the Pinatubo aerosol layer [. Mickley et al., 1997]. For example, the 

624 radiative heating at the top boundary of the aerosol layer (near 20 hPa) can lead to an 

625 enhanced, net ascent at the low latitudes. As a result, air having a relatively low mixing 

626 ratio of NOy would have been transported upward from the lower stratosphere, and the 

627 associated chemical loss of ozone due to the reactive, NOx component of the NOy would 

628 have been diminished at that time. There are similar excess ozone anomalies in the layers 

629 6U through 5U (or about 10 to 20 hPa), at least for the latitude bins of 25 S through 25N, 

630 and they lead to excessively negative linear trends when the time series analyses are 
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63 1 started at October 1991 rather than September 1 992 or January 1 993. In addition, there 

632 are greater uncertainties for the ozone signals that were attenuated by the aerosols of the 

633 lower layers of the low latitudes during this early period, so the analyses of their time 

634 series were begun at those later times, too. A further diagnosis of the true cause of the 

635 enhanced HALOE ozone in 1991-92 at 10 to 20 hPa is outside the scope of this study, but 

636 that ozone excess represents an “end point anomaly” for the determination of the trend 

637 terms (and for the coefficients of the associated 1 1 -yr terms to a lesser extent). By 

638 delaying the start date of the time series analyses by at least one year after the eruption 

639 but only for layers 6U to 4U of the lower latitudes, the effects from the volcanic aerosol 

640 layer have been reduced for the results of this present study. 

641 

642 The narrow pressure-altitude zone of rather large negative trends at 10 to 20 hPa in 

643 Figure 14 may have contributions from structures in those time series that were not 

644 accounted for. To understand that possibility, it is useful to refer to the representative 

645 zonal mean cross section of the HALOE ozone mixing ratio in Figure 6. In the low to 

646 middle stratosphere the ozone mixing ratio should be very nearly conserved during 

647 transport. Any change for its net transport at those altitudes ought to be most noticeable 

648 where the meridional gradients of the ozone mixing ratio are greatest, or in the subtropics 

649 near 10 hPa. The latitudes of 15 to 25 degrees are where the diagnosed 1 1-yr ozone 

650 responses of Figure 9 are larger. That region is also where the amplitudes of both the 

65 1 QBO and subbiennial terms are pronounced (see Figures 7 and 8). It may be that the 

652 larger negative trends and the enhanced SC-like responses in that region have a 

653 contribution from a decadal-scale, dynamical forcing. On the other hand, because the 
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654 meridional gradients for the ozone mixing ratio are nearly flat at about 20 hPa, the 

655 presence of a similar dynamical forcing would be much less apparent in the ozone at that 

656 level. There are also larger negative trends that are marginally significant at 45S and 45N 

657 and from 30 to 50 hPa. Perhaps those trends are indicating an extension of so-called 

658 “polar ozone hole” effects to lower latitudes during this 14-yr period. If that is the case, 

659 it is noted that the effects of the QBO term for the lower stratospheric ozone distributions 

660 have already been accounted for as part of the present MLR analyses. 

661 
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661 5. Summary and Conclusions 

662 This study has been based on ozone and temperature data from a single, well- 

663 characterized instrument, HALOE, and for slightly more than one complete solar cycle. 

664 An exploratory search for a solar cycle response was performed, but only for those 

665 latitudes where HALOE provided good seasonal sampling. This consideration is 

666 important for the accurate accounting of the generally, larger-amplitude seasonal terms in 

667 the time series. The HALOE ozone profiles have good S/N and a relatively high vertical 

668 resolution, both of which are important for resolving the vertical response of ozone to the 

669 relatively, small amplitude of the 1 1-yr, max minus min, uv-flux forcing. Analyses of 

670 temperature time series were restricted to the upper stratosphere and lower mesosphere, 

671 2.0 to 0.7 hPa. 

672 

673 The analyzed QBO terms for ozone have their largest amplitudes at the low latitudes and 

674 in the lower stratosphere, and their amplitudes are nearly symmetric with latitude across 

675 the two hemispheres. On the other hand, there is some asymmetry for the amplitudes of 

676 the sub-biennial terms, in accord with the small differences for its associated annual 

677 cycles in the two hemispheres. There may also be interactions between the QBO and 

678 sub-biennial terms, but their longer-period terms were not included in the analyses. 

679 

680 A simple, 1 1-yr periodic term was fit to the ozone residuals, and its phase was noted. 

681 Where the amplitude of that term was significant, its phase fit with that of the SC flux 

682 proxies, or at least within about one year of January 1991 (or 2002) — the approximate 

683 midpoint for the flux maximum of SC 22 (or SC 23). Because the effects of trend terms 
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684 can be confounded with those of the 1 1-yr terms, linear trend terms were included in the 

685 MLR analysis model. Generally, there is very good continuity for the analyzed trends 

686 with latitude and pressure-altitude, and their signs and magnitudes are reasonable for the 

687 time period of 1991-2005 from the HALOE measurements. However, it is noted that 

688 there may still be decadal-scale processes that were important, but not accounted for 

689 particularly in the upper stratosphere of the tropics and in the lower stratosphere across 

690 all latitudes. 

691 

692 There is a significant SC-like, max minus min response of order 2 to 3% in the HALOE 

693 ozone from about 2 to 5 hPa and for the latitudes of 45 S to 45N. Weaker and much less 

694 significant responses occur in the tropics near 1 hPa and from 7 to 30 hPa. There is a 

695 significant max minus min, SC-like response of order 1 K in the HALOE temperatures 

696 between 15S and 15N and from 2 to 0.7 hPa. The temperature response at mid latitudes 

697 is about half that value and much less significant. Several chemistry/climate models are 

698 yielding response profiles for both ozone and temperature that are very similar to those 

699 analyzed here [Austin et al., 2008], 

700 

701 As noted above, linear trends were diagnosed along with the other terms of the MLR 

702 analyses of the HALOE ozone, and they were near zero throughout the upper 

703 stratosphere. Negative trends of order 4 to 6%/decade were found at 10 to 20 hPa across 

704 the low to middle latitudes of both hemispheres. Weak, positive trends were found in the 

705 lower stratosphere for the subtropics of both hemispheres. Conversely, there were weak, 


33 



706 negative trends from 30 to 50 hPa in the tropics and poleward of about 25 degrees 

707 latitude of both hemispheres. 

708 

709 It is concluded that there is no real discrepancy between the SC-like response profiles 

710 from some current models and from the stratospheric ozone from the HALOE 

711 observations, at least where its SC-like terms are significant. The MLR analysis 

712 approach used herein has made it possible to separate that response from those of the 

713 expected interannual forcings and to first order from the confounding effects of any long- 

714 term trends. The precision, calibration, vertical resolution, spatial sampling , and 

715 retrieval approach of the HALOE solar occultation experiment along with the MLR 

716 analysis techniques used herein appear to be quite adequate for future determinations of 

717 the middle atmosphere responses of both ozone and temperature to the uv-flux forcing of 

718 the 1 1 -yr solar cycle. 

719 

720 
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974 Figure 1. Tangent point locations of the HALOE SR and SS measurements for 2001. 

975 

976 Figure 2. Time series of ozone (in DU) for half-Umkehr layer 8U (near 2.4 hPa) at 25N. 

977 MLR model terms are indicated at the lower left. 

978 

979 Figure 3. Time series residuals for the ozone model of Figure 2. The solid horizontal 

980 line is the least squares linear regression fit to the residuals. 

981 

982 Figure 4. Time series of temperature (in K) at 2 hPa and 25N. MLR model terms are 

983 indicated at the lower left. The solid horizontal line is the sum of its constant and linear 

984 trend terms. 

985 

986 Figure 5 . Time series of ozone for half-Umkehr layer 6L (near 13 hPa) at 25N. MLR 

987 model terms are indicated at the lower left, and the fit to the data starts in September, 

988 1992. 

989 

990 Figure 6. Zonal mean cross section of the HALOE ozone mixing ratio (in ppmv) for 

991 sunset (SS) from February 26 to April 5, 1995. Contour interval is 1 ppmv. 

992 

993 Figure 7 . Contour plot of the amplitude (in %) of the QBO term for the HALOE ozone 

994 from the MLR models. Contour interval is 0.5% from 0.0 to 5.0 and 2.0% thereafter. 

995 Darker shading is where the amplitudes exceed 2%. 

996 
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997 Figure 8. As in Figure 7, but for the amplitudes of the sub-biennial terms. 

998 

999 Figure 9. Contour plot of the max minus min 1 1-yr response (in percent) for HALOE 

1000 ozone. Contour interval is 0.5% from 0.0 to 5% and is 2% thereafter. Darker shading 

1001 denotes regions where the response exceeds a Cl of 90%; lighter shading has a Cl of 

1002 between 70 and 90%. 

1003 

1004 Figure 10. Phase variation of the 1 1-yr term of Figure 9 (in yrs from January 1991). 

1005 Contours are: dashed — negative, solid — positive, and dotted — zero, and the interval is 1 

1006 yr. The domain within the bold contours of +1 and -1 is considered as in-phase with the 

1007 SC uv-flux and that region is shaded darker. 

1008 

1009 Figure 11. Contour plot of the max minus min 1 1-yr response (in K) for HALOE 

1010 temperature. Contour interval is 0.2 K. Darker shading denotes regions where the 

1011 response exceeds a Cl of 90%; lighter shading has a Cl of between 70 and 90%. 

1012 

1013 Figure 12. Phase variation of the 1 1-yr term of Figure 1 1 (in years from January 1991). 

1014 Contours are: dashed — negative, solid — positive, and dotted — zero, and the interval is 1 

1015 yr. The domain within the bold contours of +1 and -1 is considered as in-phase with the 

1016 SC uv-flux and that region is shaded darker. 

1017 
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1018 Figure 13. Profiles of the 1 1 -yr and the SC-like, max minus min responses (in percent) 

1019 for the tropical to subtropical ozone from HALOE. The solid curve is the model result 

1 020 for 5N from Brasseur [ 1 993] . 

1021 

1022 Figure 14. Contour plot of the linear trend terms (in %/decade) from the MLR models of 

1023 HALOE ozone. Contour interval is 2%; the positive and zero contours are solid, while 

1024 the negative contours are dashed. Darker shading denotes regions where the trend has a 

1025 Cl exceeding 90%; lighter shading has a Cl of between 70 and 90%. 

1026 
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1026 

1027 

1028 Table 1 — Pressure Intervals and 

1029 Half-Umkehr Layer Designations 

1030 


Pressure Range 

Half-Umkehr 

in (hPa) 

Layer 

0.7- 1.0 

10L 

1.0- 1.4 

9U 

1.4- 2.0 

9L 

2.0- 2.8 

8U 

2.8- 4.0 

8L 

4.0- 5.6 

7U 

5.6- 7.9 

7L 

7.9-11.2 

6U 

11.2-15.8 

6L 

15.8-22.4 

5U 

22.4-31.7 

5L 

31.7-44.8 

4U 

44.8-63.3 

4L 
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1032 

1033 

1034 


1035 

1036 


Table 2 — SC-like, max minus min responses for ozone (%) and temperature (K) 
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TPA: 30.0 km 
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Latitude Progression 

Figure 1. Tangent point locations of the HALOE SR and SS measurements for 2001. 
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Figure 2. Time series of ozone (in DU) for half-Umkehr layer 8U (near 2.4 hPa) at 25N. 
MLR model terms are indicated at the lower left. 
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Figure 3. Time series residuals for the ozone model of Figure 2. The solid horizontal 
line is the least squares linear regression fit to the residuals. 
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Figure 4. Time series of temperature (in K) at 2 hPa and 25N. MLR model terms are 
indicated at the lower left. The solid horizontal line is the sum of its constant and linear 
trend terms. 
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Figure 5. Time series of ozone for half-Umkehr layer 6L (near 13 hPa) at 25N. MLR 
model terms are indicated at the lower left, and the fit to the data starts in September, 
1992. 
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HALOE ss 0 3 (ppmv) 950226 to 950405 



Figure 6. Zonal mean cross section of the HALOE ozone mixing ratio (in ppmv) for 
sunset (SS) from February 26 to April 5, 1995. Contour interval is 1 ppmv. 
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HALOE 03 QBO Amplitude (%) 



LATITUDE 


Figure 7. Contour plot of the amplitude (in %) of the QBO term for the HALOE ozone 
from the MLR models. Contour interval is 0.5% from 0.0 to 5.0 and 2.0% thereafter. 
Darker shading is where the amplitudes exceed 2%. 
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Figure 8. As in Figure 7, but for the amplitudes of the sub-biennial terms. 
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HALOE 1 1 -yr 03 Response (%) 



LATITUDE 


Figure 9. Contour plot of the max minus min 1 1-yr response (in percent) for HALOE 
ozone. Contour interval is 0.5% from 0.0 to 5% and is 2% thereafter. Darker shading 
denotes regions where the response exceeds a Cl of 90%; lighter shading has a Cl of 
between 70 and 90%. 
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Figure 10. Phase variation of the 1 1-yr term of Figure 9 (in yrs lfom January 1991). 
Contours are: dashed — negative, solid — positive, and dotted — zero, and the interval is 1 
yr. The domain within the bold contours of +1 and -1 is considered as in-phase with the 
SC uv-flux and that region is shaded darker. 
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HALOE 1 1 -yr Temp Response (K) 
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Figure 11. Contour plot of the max minus min 1 1-yr response (in K) for HALOE 
temperature. Contour interval is 0.2 K. Darker shading denotes regions where the 
response exceeds a Cl of 90%; lighter shading has a Cl of between 70 and 90%. 
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LATITUDE 


Figure 12. Phase variation of the 1 1-yr term of Figure 1 1 (in years from January 1991). 
Contours are: dashed — negative, solid — positive, and dotted — zero, and the interval is 1 
yr. The domain within the bold contours of +1 and -1 is considered as in-phase with the 
SC uv-flux and that region is shaded darker. 
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SC Ozone Response, HALOE vs. Model 



Ozone Response, Max - Min (%) 


Figure 13. Profiles of the 1 1-yr and the SC-like, max minus min responses (in percent) 
for the tropical to subtropical ozone from HALOE. The solid curve is the model result 
for 5N from Brasseur [1993]. 
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Linear Trends for HALOE Ozone (%/decade) 
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1 

2 Figure 14. Contour plot of the linear trend terms (in %/decade) lfom the MLR models of 

3 HALOE ozone. Contour interval is 2%; the positive and zero contours are solid, while 

4 the negative contours are dashed. Darker shading denotes regions where the trend has a 

5 Cl exceeding 90%; lighter shading has a Cl of between 70 and 90%. 
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